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Body  of  the  Report 

(!)•  1.  Tasks  of  the  Report 

The  Final  Report  (item  OOOAD)  must  include  a  statement  of  the  problem  studied,  and  a  summary  of  the 
most  important  results,  received  in  the  interim  reports.  The  Report  must  contain  burning  rate  and  surface 
temperature  sensitivities  of  the  studied  mixtures  (having  energetic  binders  on  the  base  of  new  oxidizer 
CL-20  (HNIW),  and  mixtures  having  HMX)  to  pressure  and  initial  temperature,  burning  rate  response 
functions  of  the  mixtures  to  linear  and  nonlinear  pulsations.  Information  about  experimental  arrangements 
and  the  obtained  temperature  profiles  and  burning  wave  parameters  at  pressure  and  sample  temperature 
variations  must  be  suggested  also.  Conclusions  about  combustion  mechanisms  of  the  studied  mixtures 
under  stable  and  pulsating  conditions  must  be  made  and  recommendations  for  future  work  must  be 
suggested. 

2.  Objects  of  Investigations. 

The  following  mixtures  of  binder/oxidizer  were  used  in  the  investigations: 

1)  GAP/HNIW,  2)  (BAMO-AMMO)/HNIW,  3)  GAP/HMX,  and  4)  (BAMO-AMMO)/HMX. 
Characteristics  of  ingredients  were  as  follows:  HNIW  (CL-20)  is  a  polycyclic  caged  nitramine,  which 
have  the  chemical  nomination  -  2,4,6,8,10,12-hexanitrohexaazaisowurtzitane;  it  is  a  powder  with  cristal 
particles  of  density  p=2,01  gr/cm^;  the  particles  have  nonregular  shape  and  sizes  from  1  up  to  150  pm; 
HMX  particles  have  sizes  less  than  200  pm;  GAP  is  an  energetic  binder,  polyglycidyl  azide,  with 
molecular  weight  about  2000  and  element  content  C36,7H5,20,6,9N4i,6;  energetic  binder  (BAMO-AMMO) 
is  a  copolymer  of  two  thermoplastic  elastomers  -  poly  bis  3,3-azidomethyloxetane  (BAMO)  and  poly  3- 
azidomethyl-3-methyloxetane  (AMMO),  in  proportion  50:50. 

Characteristics  of  the  mixtures  were  as  follows: 

1.  GAP/HNIW,  ratio  20:80,  sum  formula  023.93H2i,o5N27,96Ci7,oi;  density  p=l,78  gr/cm^  adiabatic 
(thermodynamic)  temperature  at  10  MPa  is  T„i=2927°C. 

2.  (BAMO-AMMO)/HNIW,  ratio  20:80,  sum  formula  024.2H224N27,94Ci7,2i;  p=l,76  gr/cm^  adiabatic 
temperature  at  10  MPa  is  Tad=2934°C 

3.  GAP/HMX,  ratio  20:80,  sum  formula  023.fi3H3i.7N27,67Ci6,86;  p=l,74  gr/cm^  T^2420‘’C 
at  10  MPa. 

4.  (BAMO-AMMO)/HMX,  ratio  20:80,  sum  formula  O23.9H33.05N26.95C17.21;  p=l,72  gr/cm^  T8d=2404°C  at 
10  MPa. 


3.  Microthermocouple  Techniques  and  Experimental  Arrangements. 

Temperature  profiles  of  the  combustion  waves  and  burning  surface  temperatures  were  measured  by 
microthermocouple  methods.  Ribbon  U-shaped  thermocouples  made  of  alloys  W+5%Re/W-l-20%Re  of 
2-7  pm  thick  were  imbedded  mto  the  samples.  Cylindrical  samples  of  the  mixtures  (8  mm  in  diameter  and 
3  cm  in  length)  were  burned  in  window  bomb  (volume  2000  cm^)  of  constant  pressure  in  atmosphere  of 
gas  nitrogen.  The  samples  with  2-3  thermocouples  had  been  ignited  by  electrically  heated  wire; 
thermocouple  signals  were  recorded  by  amplifier  and  oscillograph.  The  samples  were  burned  in  a  bomb 
of  constant  pressure  in  atmosphere  of  nitrogen  at  pressures  0.1  (1.0)  -10  MPa  and  at  sample  initial 
temperatures  To=  -100,  +20,  and  +100^.  Burning  rates  were  measured  by  determination  of  time  delay 
between  the  thermocouple  signals,  by  microfilming  of  sample  combustion  and  by  determination  of  time 
of  pressure  increase  during  the  sample  combustion  (pressure  change  was  measured  by  Kisstler  pressure 
detectors).  Burning  surface  temperatures  were  measured  by  thermocouples  that  are  being  pressed  to  the 
surface  during  sample  combustion  and  by  establishing  locations  of  slope  breaks  on  temperature  profiles, 
recorded  by  the  thermoeouples  (there  is  the  theory  and  experimental  substantiation  of  burning  surface 
temperature  measurements  by  determination  of  temperature  of  the  break). 


2 


The  theory  of  thermocouple  measurements  in  combustion  waves  of  solids  shows  that  aU  requirements  of 
the  theory  meet  in  thesee  investigations.  Different  types  of  metal  wires  for  thermocouples  were  used 
(We, Re  and  Pt,Rh)  to  test  the  catalytic  effect  on  thermocouples.  This  catalytic  effect  was  not  observed. 

4.  Results  of  The  Measurements. 

The  investigated  mixtures  burned  at  different  pressures  (0.5  -  10  MPa)  and  at  sample  temperatures  To 
from  -100  up  to  lOOT,  as  a  rule  relatively  steady,  without  periodically  pulsations  of  burning  rate. 
Mixture  GAP/HNIW  can  bum  also  at  0. 1  MPa. 

4.1.  Carbon  residue.  Carbon  residue  produces  on  the  burning  surface  of  mixtures  on  the  base  of  CL-20 
at  decreased  pressures.  The  residue  in  mixture  GAP/HNIW  has  increased  density  (up  to  0.05  g/cm^  at  0.1 
MPa)  and  holds  shape  of  the  sample  after  going  of  the  burning  wave  at  pressures  0.1-2  MPa.  The  density 
of  the  residue  decreases  at  pressure  increasing:  it  is  about  0.01-0.02  g/cm^  at  0.5  -2  MPa.  The  residue, 
which  holds  the  shape  of  sample,  is  a  carbon  lattice,  through  which  gaseous  combustion  products  blow 
easily.  At  pressures  2-8  MPa  the  carbon  residue  during  combustion  of  GAP/HNIW  has  forms  of  carbon 
flakes,  which  amount  decreases  with  pressure  and  at  10  MPa  the  residue  doesn’t  observed.  In  combustion 
waves  of  mixture  (BAMO-AMMO)/HNIW  the  carbon  residue  at  low  pressures  has  form  of  carbon  flakes 
and  at  elevated  pressures  the  residue  doesn’t  observed.  Carbon  residue  on  the  burning  surface  of  the 
mixtures  3,  4  (on  the  base  of  HMX)  did  not  observed.  Variation  of  To  doesn’t  influence  on  the  carbon 
residue  production. 

Results  of  experimental  measurements  of  mass  burning  rates  and  burning  surface  temperatures  in 
combustion  waves  of  the  mixtures  show  Tables  1  -  4,  which  contain  averages  values  of  these  parameters 
for  every  combustion  regime. 

4.2.  Mass  burning  rates  m.  Tables  1,  2  show  that  burning  rates  of  mixtures  on  base  of  HNIW  are  very 
high  and  that  burning  rates  of  mixture  1  are  significantly  lower  than  the  rates  of  mixture  2.  Indeed  burning 
rates  of  mixture  1  in  pressure  interval  0.5-10  MPa  comprise  0.5-4.0  g/cm^c  (To=20°C),  but  those  of 
mixture  2  -  0.8  -  6.9  g/cm^c.  It  can  be  seen  from  Tables  1  -  4  that  burning  rates  of  mixtures  on  the  base  of 
HNIW  are  significantly  higher  than  the  rates  of  mixtures  on  the  base  of  HMX.  Indeed,  burning  rates  of 
mixtures  3, 4  in  interval  0.5-10  MPa  change  from  0.1  up  to  1.25-1.45  g/cm^c  (To=20°C).  It  is  important  to 
stress,  that  influence  of  type  of  binder  on  burning  rate  of  HNIW-based  mixtures  is  significantly  higher 
than  that  on  burning  rate  of  HMX-based  mixtures. 

It  was  established  at  sample  temperature  variations,  that  burning  rates  of  mixture  2  at  negative  sample 
temperatures  have  especially  significant  changing  and  because  of  that  the  burning  rates  were  obtained  at 
-50  and  -100°C  -  see  Table  2. 

Standard  deviation  of  m  measurements  is  about  ±5  %. 

4.3.  Burning  surface  temperatures  Tg.  Tables  1  and  2  show  that  burning  surface  temperatures  Ts  of 
mixtures  1  and  2  are  significantly  lower  than  surface  temperatures  of  mixtures  3  and  4.  Indeed  the  surface 
temperatures  of  mixture  1  comprise  248-3 55°C  (0.1-10  MPa),  and  Ts  of  mixture  2  comprise  285-382°C 
(0.5-10  MPa),  wheras  surface  temperatures  of  mixture  3  comprise  330-492°C  and  330-500°C, 
correspondingly  (To~20°C).  It  can  be  seen,  that,  in  spite  of  large  burning  rates  of  mixtures  on  the  base  of 
HNIW,  the  surface  temperatures  of  these  mixtures  significantly  lower  than  values  of  Ts  of  mixtures  on  the 
base  of  HMX,  having  small  burning  rates. 

Standard  deviation  of  Ts  measurements  is  about  ±5  %. 

4.4.  Temperature  profiles.  Temperature  distributions  T(x)  along  combustion  waves  of  the  investigated 
mixtures  have  irregular  temperature  pulsations.  These  pulsations  are  significant  at  0.1-2  MPa,  possibly 
due  to  carbon  residue  creation  and  carbon  flakes  tearing  off  from  the  burning  surface.  However,  at  p>2 
MPa  the  temperature  profiles  are  relatively  smooth.  5-8  numbers  of  temperature  profiles  T(x)  were 
obtained  for  every  combustion  regime  and  then  averaged  curves  of  T(x)  were  obtained.  Figures  1-7 
present  typical  averaged  temperature  profiles  T(x),  in  solid  and  gas,  at  different  pressures  and  sample 
temperatures  for  mixture  1.  Similar  averaged  temperature  profiles  have  other  mixtures  (see  Interim  Report 
item  OOIAC).  Peculiarities  of  the  main  burning  wave  parameters  of  the  mixtures  can  be  seen  from  Tables 
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1-4,  The  figures  and  Tables  show  that  all  mixtures  have  in  the  gas  phase  a  two-zone  structure  of  the 
profiles  at  middle  pressures,  beginning  approximately  fi-om  0.5  MPa.  This  two-zone  structure  Oow- 
temperature,  or  first  zone,  and  flame  zone)  exists  up  to  about  2  MPa.  Mixture  1  can  bum  at  0. 1  MPa  and 
at  this  pressure  the  gas  phase  of  this  mixture  contains  only  one  low-temperature  zone.  At  elevated 
pressures,  beginning  from  2-6  MPa,  both  zones  merge  into  one  zone  of  the  gas  phase  for  all  mbctures  at 
all  sample  temperatures. 

Detailed  measurements  of  temperature  profiles  at  negative  sample  temperatures  were  performed. 
Mixtures  1,  3,  and  4  have  increased  irregular  temperature  pulsations  at  To=  -lOOT.  Mixture  2  showed 
more  re^lar  temperature  pulsations  at  To=  -100°C  at  p<5  MPa.  It  gives  some  base  for  suspicion  that 
combustion  of  this  mixture  under  indicated  condition  has  unstable  character. 

5.  The  Experimental  Data  Processing. 

The  processing  of  the  obtained  temperature  profiles  allows  a  set  of  burning  wave  parameters  to  be 
detennined.  Tables  1-4  present  the  dependencies  of  the  parameters  on  pressure  and  To.  Table  5  shows 
estimated  due  to  gas  product  content  the  dependencies  of  coefficients  of  heat  conduction  Xg  and  specific 
heat  Cp  of  the  gas  phase  on  temperature. 

5.1.  Temperature  and  geometric  characteristics  of  temperature  distributions  in  gas. 

Values  of  mean  temperature  of  the  first  zone,  Ti,  for  mixture  1  are  equal  to  800-1 200'’C  and  for  mixture  2 
-  1000-1 150°C  (To=20‘’C).  Temperatures  of  Ti  for  mixture  3  comprise  980-1050°C  and  for  mixture  4  - 
900-1020'’C.  It  was  estblished  that  sample  temperature  variations  change  T|  approximately  by  similar 
way.  Standard  deviation  of  T|  measurements  is  about  ±10  %, 

Flame  temperatures  Tf  for  mixture  1  are  equal  to  2500-2930°C  (To=20‘’C)  and  for  mixture  2  -  2200- 
2930'’C  (To=20"C).  Values  of  Tf  for  mixture  3  comprise  1 900-2420®C  and  for  mixture  4  -  2000-2400°C.  It 
can  be  seen  that  flame  temperature  of  mixtures  on  base  of  HNIW  is  much  higher  than  that  of  mixtures  on 
base  of  HMX.  Maximal  flame  temperatures  Tf  were  obtained  at  5-8  MPa.  It  was  estblished  that  sample 
temperature  variations  change  Tf  approximately  by  similar  way. 

Distances  Li  between  the  burning  surface  and  the  beginning  of  the  flame  zone  (L|  is  the  length  of  the 
first  flame)  decrease  when  pressure  increases.  Values  of  Li  for  mixture  1  are  equal  to  2-  1  mm  (0.1-2 
MPa),  for  mixture  2  -  to  1-0.25  mm  (0.5-5  MPa),  for  mixture  3  -  to  0.6-0.5  mm  (0.5-1  MPa),  and  for 
mixture  4  -  to  0.5-0. 4  mm  (0.5-1  MPa)  -  all  for  To=20°C.  It  can  be  seen  that  mixtures  on  base  of  HNIW 
have  greater  length  of  the  first  flame  than  mixtures  on  base  of  HMX. 

Distances  L  between  the  burning  surface  and  the  flame  (up  to  ~0.99  Tf)  also  decrease  when  pressure 
increases.  In  fact,  L  is  the  length  of  the  gas  phase  reaction  zone.  Values  of  L  for  mixture  1  are  equal  to  4- 
1.8  mm  (0.5-10  MPa),  for  mixture  2  -  to  5-2  mm  (0.5-10  MPa),  for  mixture  3  -  to  2-1.3  mm  (0.5-10  MPa), 
and  for  mixture  4  -  to  2. 2-0. 9  mm  (0.5-10  MPa)  -  all  for  To=20‘’C.  It  can  be  seen  that  mixtures  on  base  of 
HNIW  have  larger  length  of  the  gas  phase  reaction  zone  than  mixtures  on  base  of  HMX. 

Standard  deviation  of  Li  andL  measurements  are  ±(15-30)%. 

5.2.  Coefficients  of  heat  diffusivity  and  heat  conduction  in  solid  heat  layer. 

Coefficients  of  solid  heat  diffusivity  were  obtained  by  formula  x—lfb,  where  /  is  thickness  of  heat 
layer  of  the  condensed  phase  (see  next  section).  Tables  1-4  show  that  values  of  x  increase  when  pressure 
increases.  The  values  of  %  for  mixtures  1 , 2  increase  from  2- 1 0'^  cmVs  up  to  6- 1 0’^  cmVs.  For  mixtures  3, 4 
the  values  of  %  comprise  (1.0-3.8)10‘^  cmVs.  These  values  are  typical  ones  for  %  of  solid  heat  layer  of 
many  propellant  mbctures. 

The  coefficients  of  solid  heat  conduction  were  obtained  by  formula  X.s=x-c  p.  It  was  assumed  that 
solid  specific  heat  is  equal  to  c=0.35  cal/g  K  for  mixtures  1  and  2  and  c=0.4  cal/g  K  for  mixtures  3, 4. 
Tables  1,  2  show  that  Xs  increases  with  pressure,  and  increases  when  mean  solid  heat  layer  temperature 
increases.  Values  of  X,  comprise  (10-40)10’^  cal/cm  s  K  for  mixtures  1,  2  and  (7-26)10"*  cal/cm  s  K  for 
mixtures 
3,4. 
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5,5.  Characteristic  sizes  in  solid  and  the  gas  phase. 

Characteristic  thickness  /  and  Im  in  solid,  presented  on  Tables  1-4,  were  obtained  by  temperature 
distributions  in  the  condensed  phase.  Experimentally  obtained  thickness  of  heat  layer  in  sohd  /  is  a 
distance  between  the  burning  surface  and  a  section  in  sohd  with  temperature  T*=(Ts-To)/e  +  To 
(e  is  the  base  of  In).  Thickness  of  L  of  the  partly  melt  layer  (with  melt  oxidizer)  in  the  condensed  phase  is 
a  distance  between  the  burning  surface  and  a  section  with  oxidizer  melt  temperature  Tm.  It  was  assumed 
that  the  melt  temperature  for  CL-20  is  equal  to  Tin=190°C  and  for  HMX  -  Tm=280°C.  Tables  show  that 
thickness  of  1  and  Im  decrease  when  pressure  increases.  Values  of  L  are  always  less  than  values  of  /, 
especially  for  mixtures  3,4. 

Conductive  sizes  ^  of  the  gas  phase  were  obtained  by  formula:  0=Xg/mCp;  where  Xg  is  coefficient  of  heat 
conduction  in  gas  and  Cp  is  coefficient  of  specific  heat  of  the  gas  phase.  Values  and  Cp  were  estimated 
by  content  of  combustion  products.  The  obtained  dependencies  of  Xg(T)  and  Cp(T)  are  presented  on  Table 
5.  Values  of  Xg,  which  were  used  for  S  calculations,  were  taken  at  temperatures  Ti  for  two-zone  structures 
and  at  Tf/2  for  one-zohe  structures.  Tables  1-4  show  that  values  of  0  quickly  decrease  when  pressure 
increases.  Values  of  0  comprise  18-1.8  pm  for  mixtures  1, 2  and  51-4  pm  for  mixtures  3, 4.  It  can  be  seen 
that  always  the  following  inequalities  take  place:  d«Li  and  0«L. 

Relative  thickness  of  the  gas-phase  reaction  layer  were  estimated  by  formula  Q-Li/S  MPa  for  two-zone 
structures,  or  by  formula  Q=L/0  for  one-zone  structures.  It  can  be  seen  that  inequality  of  Q»1  always 
takes  place.  It  implies  that  a  vvdde  reaction  zone  in  the  gas  phase  is  observed  for  all  combustion  regimes  of 
combustion  of  mixtures  1-4. 

5.4.  Heat  manifestations:  heat  release  in  solid  and  heat feedback  from  gas  into  solid. 

The  temperature  profiles  and  burning  surface  temperatures  were  used  for  obtaining  heat  parameters  of  the 
burning  waves.  Heat  flux  fi-om  gas  to  solid  q-m  by  heat  conduction  was  estimated  by  the  foUowing 
formula: 

q.m=-Xg(T)-(dT/dx)o;  (5.1) 

Here:  (dT/dx)o=cp  is  the  temperature  gradient  in  gas  close  to  the  burning  surface.  Values  of  cp  were 
determined  by  slopes  of  the  temperature  distributions  T(x)  in  gas  close  to  the  surface. 

Heat  feedback  q  from  gas  into  solid  by  heat  conduction  was  estimated  by  the  following  formula: 

q=-VT)-cp/m;  (5.2) 

Values  of  Xg  were  taken  at  temperatures  T=Ts+100°C. 

Formula  for  heat  release  Q  in  reaction  layer  of  the  condensed  phase  (or  on  the  burning  surface) 
have  the  following  forms: 

for  mixtures  1, 2  Q=c  (Ts-To)-q-qr ;  (5.3) 

and 

for  mixtures  3, 4  Q=c  (Ts  -To)-q-qr+qm;  (5 .4) 

where  q^  is  the  heat  of  HMX  melting  in  mixtures  3,  4;  heat  of  HNIW’s  melting  are  not  known  and 
because  of  that  formula  for  mixtures  1, 2  don’t  have  member  of  qm.  Heat  of  HMX’s  melting  is  equal  to 
24  cal/g  (per  gram  of  HMX)  and  because  of  that  qm=19.2  cal/g  (per  gram  of  mixtures  3,  4).  Radiant  heat 
feedback  qr  from  flame  into  sohd  was  estimated  for  all  mixtures  by  using  of  content  of  combustion 
products  and  real  flame  thickness  of  the  burning  samples  in  the  bomb  of  constant  pressure.  The  heat 
rachations  of  CO2,  H2O,  H2,  etc  were  taken  into  considerations.  It  was  assumed  that  coefficient  of  heat 
radiation  absorption  of  the  burning  surface  was  equal  to  0.5.  Tables  1-4  show  the  obtained  heat 
parameters. 

It  can  be  seen  from  Tables  1-4  that  all  investigated  mktures  have  significant  positive  heat  release  in  solid 
Q.  Parameter  Q  comprise  at  pressures  0.5-10  MPa  the  following  values:  60-140  cal/g  for  mixture  1,  60- 
140  cal/g  for  mixture  2, 78-180  cal/g  for  mixture  3,  and  96-185  cal/g  for  mixture  4.  Parameter  Q  increases 
always  with  pressure  and  decreases  with  sample  temperature  increasing.  It  can  be  seen  that,  as  a  rule,  the 
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values  of  Q  for  HMX-based  mixtures  higher  than  for  HNlW-based  mixtures.  Temperature  gradient  <p  in 
gas  close  to  burning  surface  has  significant  values  {(2-15)-10'^,  K/cm}  and  always  increases  with  pressure 
and  decreases  with  sample  temperature  increasing.  On  the  contrary,  heat  feedback  q  fi-om  gas  into  solid 
always  decreases  when  pressure  increases.  Values  of  q  are  comparable  with  values  of  Q  only  at  low 
pressure  (0.5  MPa)  and  only  for  HMX-based  mixtures.  Under  other  conditions  values  of  q  are 
significantly  smaller  than  values  of  Q,  especially  at  elevated  pressures.  Indeed,  Tables  1,  2  show  that  q 
comprises  in  the  studied  pressure  intervals  the  following  values;  1 1-6.3  cal/g  for  mixture  1, 9-4.3  cal/g  for 
mixture  2, 64-25  cal/g  for  mixture  3,  and  46-22  cal/g  for  mixture  4.  At  last,  the  values  of  q,  are  small;  they 
comprise:  4-8  cal/g  for  mixture  1,  2. 5-4. 6  cal/g  for  mixture  2, 0.7-3. 7  cal/g  for  mixture  3,  and  0.6-4  cal/g 
for  mixture  4.  It  can  be  seen,  that  total  heat  feedback  q+qr  fi'om  gas  to  solid  surface  less  than  heat  release 
in  sohd,  i.e.  inequality  Q>q+qr  always  take  place. 

5.5.  Heat  release  rate  in  gas  close  to  the  burning  surface. 

Heat  release  rates  Wo  in  gas  close  to  the  surface  were  estimated  by  using  of  the  obtained  data.  Wo  shows, 
obviously,  the  intensity  of  chemical  reactions  in  gas.  Values  of  Wo  were  obtained  by  the  following 
expression; 

Wo=Cp-m-(p;  (5.5) 

This  expression  was  obtained  fi'om  equation  of  heat  conduction  for  stable  propagating  combustion  waves 
for  the  case  of  significant  positive  heat  release  in  solid  and  wide  reaction  zone  in  gas.  These  experimental 
facts  were  established  in  the  previous  sections  and  it  implies  that  the  indicated  expression  for  Wo  is  valid. 
Values  of  Cp  here  were  taken  at  temperatures  T^T^-f  100°C.  Tables  1,  2  show  that  values  of  Wo  quickly 
increase  with  pressure.  The  heat  release  rate  Wo  comprise  at  0.5-10  MPa  the  following  values;  5-177 
kcal/cm^s  for  mixture  1,  10.8-355  kcal/cm^s  for  mixture  2,  1.4-64  kcal/cm^s  for  mixture  3,  and  0.9-73 
kcal/cm  s  for  mixture  4  (at  To=20^C).  It  can  be  seen  that  HNIW-based  mixtures  have  more  intensive 
chemical  reactions  in  gas  near  the  burning  surface  than  mixtures  on  the  base  of  HMX. 

5. 6.  Burning  Rate  Control  Regions, 

The  results  of  measurements  of  parameters  q+qr  and  Q  and  the  investigations  of  the  gas  phase 
zone  behavior  show  that  burning  rate  of  the  mixtures  are  caused  by  the  heat  release  in  solid  (or  on  solid 
surface)  and,  in  a  smaller  degree,  by  the  heat  feedback  from  the  gas  layer  which  is  close  to  the  burning 
surface.  It  was  mentioned  above  also  that  the  values  of  q+qr  decreases  when  the  burning  rate  increases. 
These  facts  show  that  heat  release  in  solid  (or  on  solid  surface)  is  a  main  factor  of  creating  of  the  burning 
rate  of  the  mixtures.  Heat  feedback  from  gas  to  solid  is  only  small  additional  factor  of  the  burning  rate 
creating  and  influence  of  this  factor  on  values  of  m  decreases  when  pressure  increases  (at  pressures  p^IO 
Nffa).  Then,  the  obtained  data,  presented  in  Tables  1-4,  allow  an  important  conclusion  to  be  made:  the 
high  temperature  region  of  the  burning  waves  cannot  affect  the  burning  rate.  Indeed,  the  influence  of  the 
heat  release  in  the  gas  phase  on  the  burning  rate  can  be  estimated  by  the  following  formula  obtained  as  a 
solution  of  the  heat  conduction  equation: 

00 

m  q=J  <l)(x)  exp(-x/9)  dx;  (5.6) 

0 

This  formula  shows  that  influence  of  heat  release  in  gas  (here  <I>(x)  is  distribution  of  heat  release  rate  in 
the  gas  phase)  on  m  decreases  veiy  quickly  when  x>9.  Tables  1-4  show  that  values  9  are  small  and  they 
are  sigmficantly  smaller  than  gas  zone  sizes  L  and  Li:  see  large  values  of  Cl  for  all  combustion  regimes. 
That  implies  that  high  temperature  region  of  the  burning  waves  cannot  affect  the  burning  rate.  Thus  the 
obtained  results  show  that  the  burning  rate  control  regions  in  the  combustion  waves  for  all  regimes  of  the 
mixtures  combustion,  under  investigated  conditions,  are  the  regions  of  heat  release  in  solid  just  under  the 
burning  surface  (or  immediately  on  the  burning  surface)  and  thin  low-temperature  gas  layers  close  to  the 
burning  surface  (at  x«9).  High  temperature  gas  regions  cannot  influence  Ae  burning  rate  in  fact  because 
of  a  very  large  heat  resistance  of  the  gas  phase  of  the  mixtures. 

5. 7.  Macrokinetics  of  Solid  Gasification. 

The  follo\<ing  equation  connects  burning  rate  of  solid  with  bum-surface  temperature  and 
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macrokinetic  characteristics  of  solid  gasification: 


m^=Xsp/Q^-S«Ts^/E-koQ*-l/N-exp(-EmTs);  (5.7) 

where  N=l/ris+(l-Tis)/Ti.-ln(l-Tis)-(q/Q)-[ln(l-Ti,)]/ris;  ris=Q/Q*;  Q*  is  the  maximum  of 
the  heat  release  in  solid,  E  is  activation  energy  of  limiting  stage  of  the  gasification  process,  Ts  in  K.  For 
studied  combustion  waves  Nsl.  The  expression  was  obtained  by  solution  of  system  of  two  equations  for 
the  steady  propagated  burning  wave  (the  heat  conduction  equation  for  soUd  phase  and  the  equation  of 
diffusion  of  reagents  were  solved).  The  burning  wave  propagates  due  to  heat  release  in  solid  Q  and  heat 
feedback  fi-om  gas  to  sohd  q.  Function  of  the  volumetric  heat  release  rate  Os  in  sohd  reaction  layer  was 
assumed  as  follows: 

Os(ti,  T)=Q*ko-p-(l-Ti)-exp(-E/5lT);  (5.8) 

Here  t]  is  the  reaction  completeness  (11=113  on  the  burning  surface),  ko  is  the  pre-exponent  multiplier.  The 
member  exp(-E/5RTs)  in  the  equation  above  plays  the  most  important  role  in  this  formula.  Because  of  that 
the  connection  m  and  Ts  can  be  presented  by  the  following  simpUfied  expression: 

m=A-exp(-E/25RTs);  (5.9) 

This  expression  has  been  designated  “gasification  law”.  The  results  of  measurement,  presented  in  Tables 
1,  2,  allow  the  following  gasification  law  for  mixtures  1  and  2  (i.e.  for  mixtures,  based  on  HNIW)  to  be 
obtained: 

m=1.810®exp(-32800/23lTs);  (m in g/cm\ Ts  in K)  (5.10) 

Possibly,  value  E=32800  kcal/mole  is  an  activation  energy  of  the  limiting  stage  of  the  process  of  oxidizer 
CL-20  gasification  in  the  reaction  layer  of  mixtures  1  and  2. 

The  results  of  measurement,  presented  in  Tables  3, 4,  allow  the  following  gasification  laws  for  mixtures  3 
and  4  to  be  obtained: 

m=l.M0"-exp(-28000/29lTs);  (m  in  g/cm\  T  in  K);  (5.11) 

(vahdatTs<500“C), 

The  obtained  value  of  activation  energy  of  the  gasification  process  is  very  close  to  that  of  cychc 
nitramines  (when  Ts<500°C).  Possibly,  gasification  of  oxidizer  HMX  m  the  reaction  layer  of  burning 
mixtures  3  and  4  is  the  limiting  stage  of  the  process  of  the  mixtures  gasification. 

6.  Pressure  and  Temperature  Sensitivities  of  Burning  Rates 
and  Surface  Temperatures  and  Criteria  of  Stable  Combustion. 

6.1.  The  Pressure  and  Temperature  Sensitivities. 

The  following  pressure  and  temperature  sensitivities  for  burning  rates  and  burning  surface  temperatures 
were  foimd  from  the  data  presented  on  Tables  1-4: 

P=(51nm/5To)p-coiist,  -  temperature  sensitivity  of  mass  burning  rate; 

r=(5Ts/5To)p.oonst,  -  temperature  sensitivity  of  burning  surface  temperature; 

v=(01nin/51np)7'o-oonst,  -  pressure  sensitivity  of  mass  burning  rate; 

p.=(Ts-To)‘’-(5Ts/91np)ro-oonst,  -  pressure  sensitivity  of  burning  surface  temperature; 

The  obtained  results  for  the  investigated  mixtures  are  shown  on  Tables  6  for  To=20‘'C.  The  Table  also 
contains  parameter  k,  which  was  found  by  formula  k=P  (Ts-To). 

It  can  be  seen  that  bum-rate  pressure  sensitivities  v  for  mixture  1  increase  with  pressure  :  parameter  v 
is  equal  to  0.37-0.6  at  pressures  0.1-2  MPa,  and  v=0.8-0.93  at  pressures  5-10  MPa.  Mixture  2  has 
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increasing  values  of  v  at  0.5-2  MPa  (v=0.55-0.82)  and  decreasing  v  atp>2  MPa  (v=0.73-0.68).  Parameter 

V  decreases  with  pressure  for  mixtures  3  and  4  from  v=0.83-0.9  at  1  MPa  to  v=0 .76-0.79  at  10  MPa. 

It  can  be  seen,  fliat  bum-rate  pressure  sensitivities  v  for  HMX  -based  mixtures,  as  a  mle,  are  higher  than 

V  for  HNIW-based  mixtures  (an  exclusion  for  values  of  v  for  mixture  1  at  p>5  MPa,  which  are 
extraordinary  high). 

Standard  deviation  of  values  v  determinations  is  ±  (1-3%). 

Temperature  sensitivity  of  mass  burning  rate  p  has  a  typical  dependence  on  pressure:  a  maximum  at 
0.5-2  MPa  and  decreasing  with  pressure  increasing.  Values  of  p  are  close  each  other  for  mixtures  1, 3, 4. 
Indeed,  p=0.23-0.28  %/K  at  I  MPa  and  decreases  to  0.12-0.17  %/K  at  10  MPa  for  the  mixtures.  Values  of 
P  for  mixture  2  are  much  higher  than  that  of  other  mixtures.  Indeed,  P=0.3  %/K  at  1  MPa,  0.36  at  2  MPa 
and  0.28  at  10  MPa. 

Standard  deviation  of  values  p  determinations  is  ±  (5-6%). 

Parameters  of  k  for  mixtures  1,  3, 4  are  less,  than  1  at  all  pressures,  only  for  mixture  2  at  p  >  2  k  >  1. 
Very  low  values  k  has  mixture  1 :  k  <  0.6  at  1  - 1 0  MPa. 

Standard  deviation  of  values  k  determinations  is  ±  10%. 

Temperature  sensitivity  of  surface  temperature  r  for  HMX-based  mixtures  1  has  decreasing  dependence 
on  pressure  and  comprise  r=0.17-0.19  at  1  MPa  and  decreases  to  1=0.09-0.067  at  10  MPa.  Mixture 
GAP/HNIW  has  lower  values  of  r:  this  parameter  is  almost  constant  in  pressure  interval  0. 1-10  MPa  and 
comprises  0.06-0.1.  Mixture  (AMMO-BAMO)/HNrW  has  increasing  dependence  of  r  on  pressure  and 
this  parameter  is  equal  to  0. 13  at  0.5  MPa  and  increases  to  0.2  at  10  MPa.  Mixture  2  has  maximal  values 
of  r  at  p  >  2  MPa. 

Standard  deviation  of  values  r  determinations  is  ±  1 0%. 

Pressure  sensitivities  of  burning  surface  temperature  |x  for  HMX-based  mixtures  are  close  each  other 
and  comprise  0.16-0.17  at  1  MPa  and  decrease  to  0.09  at  10  MPa.  HNIW-based  mixtures  have  other 
dependencies  for  p  on  pressure:  parameter  p  of  mixture  2  practically  independent  on  pressure  and 
comprises  0.1-0.12,  and  p  of  mixture  1  increases  from  0.06  at  0.1  MPa  to  0.14  at  10  MPa. 

Standard  deviation  of  values  p  determinations  is  ±  (1-3%). 

The  obtained  sensitivities  are  the  most  important  characteristics  of  combustion  processes  of  solid 
mixtures.  These  characteristics  are  used  for  characterizations  of  solid  mixtures  at  stable  and  at  pulsating 
combustion.  The  obtained  above  sensitivities  allow  ones  to  show  that  well-known  Jacobian  5  =  v-r-p-k  is 
very  close  to  zero.  Indeed,  accuracy  of  5  obtaining  is  about  ±0.005,  and  it  is  inside  of  experimental 
accuracy  of  the  determination.  It  implies  that  equation  (5.9)  has  experimental  validation. 

Estimation  of  criteria  of  stable  combustion  (see  next  section)  is  the  second  example  of  usage  of 
the  sensitivities  for  determination  of  region  of  stable  combustion. 

6.2.  Criteria  of  Stable  Combustion. 

The  simplest  criterion  of  stable  combustion  is  Zeldovich’s  criterion  k=p  (Ts-To).  It  is  used,  when  surface 
temperature  is  changing  no  significantly  at  pressure  and  sample  temperature  variations.  The  criterion 
predicts  the  stable  combustion  at  k  <  1.  However,  practically  always,  the  significant  changing  of  the 
surface  temperature  takes  place  (see  Tables  1-4).  In  this  case  at  k  >  1  (at  k  <  1  Zel’dovich’s  and 
Novozhilov’s  criteria  are  identical)  Novozilov’s  criterion  k*,  which  takes  into  consideration  T*  changing, 
must  be  used.  The  criterion  predicts  stable  combustion  at  k*  <  1.  The  formula  for  the  criterion  is  as 
follows: 

k»=(k-l)V(k+l)r;  (6.1) 

Calculated  criterion  k  for  the  investigated  mixtures  shows  Table  6.  It  can  be  seen  that  mixtures  1, 3,  and  4 
have  k  <  1.  It  shows  that  these  mixtures  have  stable  combustion  under  indicated  conditions.  Further 
estimations  showed  that  even  at  T„=  -100®C  mixtures  1, 3,  and  4  have  stable  combustion  because  for  these 
mixtures  k*  <  1  at  all  investigated  sample  temperatures  (at  To  =  -1 OOX)  mixtures  1, 3,  and  4  have 
k  >  1  at  low  pressures). 

Table  6  shows  that  at  To=20“C  mixture  (AMMO-BAMO)/HNIW  has  k  >  1  at  pressures  2-10  MPa. 
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But  estimations  of  k*  values  show  that  this  mixture  has  also  stable  combustion  even  at  To  >  -50°C  -  see 
Table  7.  However,  at  To  =  -100°C  and  p  <  5  MPa  for  mixture  2  inequality  k*  >  1  takes  place.  It  impUes 
that  at  To  =  -100“C  and  p  <  5  MPa  mixture  (AMMO-BAMO)/HNIW  bums  unstable.  This  conclusion 
confirms  results  obtained  for  temperature  profile  measurements  -  see  above,  section  4.4. 


7.  Second  Step  of  Data  Processing  ; 

Combustion  Characterizations  at  Pulsating  Pressure 

7.i.  Pressure-Driven  Burn-Rate  Response  Functions  at  Linear  Approximation. 

The  obtained  results,  which  were  received  at  constant,  nonpulsating  pressures,  allow  us  to  study  burning 
rate  behavior  of  the  investigated  mixtures  under  conditions  of  oscillatoiy  pressure.  It  can  be  made  due  to 
obtained  above  sensitivities  form  and  T*.  This  paragraph  is  devoted  to  this  study. 

Let  us  assume  that  external  pressure  changes  by  harmonic  way  and  is  equal  to  : 

p=Po+prCOSC3f,  (7.1) 

here  po'is  a  constant  mean  pressure  level  and  pi  is  an  oscillating  pressure  (pi«Po),  t  is  current  time,  and  tn 
is  cyclic  firequency  of  oscillation  in  1/s. 

Then  mass  burning  rate  in  linear  approximation  is  equal  to: 

m=mo+miCOs(Bjt+(^) ;  (7.2) 

where  m  is  current  burning  rate,  nio  is  stationary  burning  rate  at  pressure  po ,  mi  is  a  pulsating  binning  rate 
(mi«m),  and  ifr  is  phase  shift. 

Then  dimensionless  burning  rate  has  the  following  form; 

F=l+Fiocos(sr+  (^),  (7.3) 

here  F=m/mo ,  Fio=mi/m<,. 

It  is  convenient  to  analyze  burning  rate  response  in  complex  form.  In  this  case  complex  dimensionless 
burning  rate  is  equal  to  F=l+Fio-e  '®'  and  complex  dimensionless  pressure  is  equal  to  where 

Tli=Pi/po,  and  CO  is  dimensionless  frequency. 

The  following  expression  for  burning  rate  response  fimction  to  harmonic  oscillatory  pressure  at  linear 
approximation  takes  place: 

U=  [v  +  (v  r  -  p  k)  (z  - 1)]/[1  +  r  (z  - 1)  -  k  (z  -  l)/z];  (7.4) 

Here  z=(l + Vl  +  4/o  )/2;  U=Vio/?7i,  Vio  is  dimensionless  burning  rate  complex  amphtude,  and  co  is 
dimensionless  frequency,  which  is  equal  to  the  cychc  frequency  in  1/s,  multiphed  by  the  time  of  solid 
heat  layer  relaxation  x/fb^  Values  of  v,  p,  k  and  r  here  are  based  on  the  mean  pressure  po  and  because  of 
that  data  from  Table  6  can  be  used  for  estimations  of  complex  fimctions  of  U.  Thus,  indeed,  we  can  study 
the  burning  rate  behavior  of  the  investigated  mixtures  under  conditions  of  pulsating  pressure.  The 
equation  (7.4)  is  based  on  the  supposition  that  relaxation  time  of  sohd  heat  layer  is  much  longer  that  those 
for  reaction  layer  of  sohd  and  the  leading  part  of  the  gas  phase.  Estimations  show  that  for  the  investigated 
mixtures  this  supposition  is  valid:  values  of  %lx^  are  at  least  10-20  times  higher  than  other  relaxation 
times.  Because  of  that,  all  estimations  of  the  response  functions  were  made  for  dimensionless  frequency 
interval  co=0-20. 

Typical  dimensionless  amplitudes  Re{U}of  burning  rate  response  fimctions  and  phase  of  the  response, 
Im{U},  in  radians,  as  dependencies  on  co  are  shown  on  Figures  8-14  for  mixture  1.  Similar  functions 
were  obtained  for  other  mixtures.  Pecuharities  of  characteristics  of  the  functions  show  Table  8  (see 
below).  It  can  be  seen  that  the  real  part  of  U(<a)  may  reach  values,  which  are  much  greater  than  its  steady 
state  value  Re{U(0)}=v.  It  is  a  consequence  of  high  quality  (in  filter  theory  terminology)  of  our 
oscillatory  systems.  The  maximum  of  Re{U}  is  at  a  frequency,  which  are  close  to  the  natural  (resonance) 
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fiequency  cv=  4k  /r.  In  this  frequency  area  the  denominator  of  equation  (7.4)  is  less  than  unity.  Thus, 
any  small  perturbations  of  the  pulsating  system  can  significant  influence  the  response  frmction.  It  is 
natural,  that  function  Im(U(o)},  presenting  the  phase  shift  of  the  response  function,  changes  the  sign  at 
the  resonance  fixquency. 

Table  8  shows  results  of  calculations  of  characteristic  points  of  Re{U}  and  Im{U}  for  the  studied 
mixtures.  Relative  standard  deviations  of  values  Re{U}m  obtaining  are  ±  (10-30%). 

The  used  denominations  in  Table  8  are  as  follows: 

Re(U}ni  is  maximum  of  Re{U(«a=a)„)}  in  interval  o  =0-20;  -  dimensionless  frequency,  corresponding 

to  this  maximum  (the  resonance  frequency);  Re{U}2o  =  Re{U(20)};  Im{U„4  -  maximum  of  Im{U}  in 
interval  a  =0-c^,  i.e.  in  positive  region  of  values  of  ImU;  0  m  ~  dimensionless  frequency,  corresponding 
to  this  maximum;  Im{Um.}  -  minimum  of  ImU  (if  this  minimum  takes  place)  in  negative  region  of  ImU; 

Table  8  shows  that  response  functions  Re{U(fii)}  for  HNIW-based  and  for  HMX-based  mixtures 
weakly  depend  on  0,  especially  at  elevated  p.  The  result  shows  that  resonance  phenomena  in  the  studied 
burning  mixtures  are  significantly  depressed  due  to  viscous  behavior  of  the  oscillations.  In  fact,  the 
pulsating  combustion  of  the  mixtures  is  an  example  of  a  relatively  smoothed  resonance  phenomenon. 

Table  8  shows,  that  changing  of  oxidizer  HMX  to  oxidizer  HNTW  in  mixture  with  binder  GAP,  improves 
pulsating  behavior  of  burning  rate.  Indeed,  mixture  GAP/HMX  has  Re{U}„=I.28-1.56,  and  mixture 
GAP/HNIW  has  Re{U}n,=0.8-1.3.  The  obtained  result  is  very  important,  because  previous  investigations 
showed  that  increasing  of  energy  of  mixtures  leads,  as  a  rule,  to  deterioration  of  pulsating  behavior  of 
burning  rate.  But  GAP/HNIW  has  higher  energy,  than  GAP/HMX.  So,  HNIW  in  mixture  GAP/HNIW 
demonstrates  exceptions  to  this  rule.  It  is  an  advantage  of  the  new  oxidizer  in  mixture  GAP/HNIW. 

On  the  contrary,  mixture  (AMMO-BAMO)/HNIW  has  increased,  in  comparison  with  mixture 
(AMMO-BAMO)/HMX,  values  of  Re{U}m.  Indeed,  mixture  2  has  Re{U}„=I.2-2,  i.e.  high  response.  But 
mixture  4  has  Re{U}m=I.2-1.6.  Thus,  binder  (AMMO-BAMO)  doesn’t  suitable  for  using  HNIW  as  an 
oxidizer  with  binder  AMMO-BAMO,  at  any  rate  without  catalytic  additives. 

So,  mixture  1  is  the  best  mixture,  having  smallest  linear  response  of  burning  rate  on  acoustic  pressure 
pulsations;  mixture  4  occupies  the  second  place,  on  the  third  place  stands  mixture  3,  and  the  worst  one  is 
mixture  2.  This  situation  for  values  Re{U}m  of  the  studied  mixtures  in  linear  approximation  can  be 
presented  as  following  «  sequence  of  preferences  »  ; 

[GAP/HNIW]  <  [(AMMO-BAMO)/HMX]  <  [GAP/HMX]  <  [(AMMO-BAMO)/HNIW]  -  LIN.  (k) 

On  the  left  side  here  stands  the  minimal  response  and  on  the  right  side  -  the  maximal  one.  Mark  ‘LIN’ 
show,  fliat  it  is  the  linear  case.  It  is  interesting  to  note,  that  difference  [Re{U}m  -  u]  also  corresponds  to 
inequality  (k). 


Pressure-Driven  Burn-Rate  Response  Functions  at  Nonlinear  Approximation. 

The  Second  Order  Interactions. 

This  section  introduces  the  defimtion  of  response  function  of  the  higher  orders.  In  a  common  case  the 
acoustic  field  contains  several  harmonics  (modes).  No  linearity  of  combustion  process  can  give  the 
oscillation  of  the  burning  rate  with  different  frequencies  by  interaction  of  two  (or  several)  other  modes. 
The  simplest,  and  more  significant,  are  the  following  cases; 

1)  Interaction  of  the  first  and  the  second  harmonics  gives  the  second  order  correction  to  linear 
response  fimetion  for  the  first  and  third  harmonics  (because  the  product  coscox  x  cos2o)x  contains 
oscillations  with  co  and  3<o  frequencies), 

2)  Self-interaction  of  the  first  mode  gives  the  second  mode  (because  the  product  cosox  x  cosox 
contains  oscillations  with  2<a  frequency). 

Such  nonlinear  interaction  may  produce  a  significant  effect  in  the  case  when  one  of  interacting  harmonics 
is  close  to  the  natural  frequency  oh. 

Values  of  v,  p,  k  and  r  here  are  based  also  on  the  mean  pressure  po,  as  in  the  linear  case. 

In  the  next  sections  for  nonlinear  approximations  the  italic  letters  were  used  for  the  same  parameters:  v, 
nk,r,po,p,p\. 
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7.2.1.  The  second  order  interaction  of  the  first  and  the  second  harmonics. 

It  was  mentioned  above,  that  acoustic  field  contain  in  a  common  case  several  harmonics.  For  the  second 
order  interactions  of  the  harmonics  we  must  use  the  following  expressions: 

p=Pa  +  jjicosor + p2COs{l(or+  (7.5) 

m=mo  +  miCOs((ar)-\i/mi)  +  m2cos(20rf\(/m2);  (7 .6) 

After  presenting  of  the  expressions  for  dimensionless  pressure  and  burning  rate  in  complex  forms  and 
solving  of  the  no  stationary  heat-conduction  equation  in  the  condensed  phase,  the  following  solution  can 
be  obtained  for  linear  response  function  plus  nonlinear  correction  of  fee  fimction  by  the  second  order 
interaction  of  the  first  and  second  modes: 


V|  —  U  +  Un-rintjfy}! ; 


(7.7) 


where  :  U  -  linear  response  fimction,  see  eq.(7.4),  and  Unpntjffji  is  the  member  of  the  nonlinear 
correction.  Here : 


Un  = 


gi 


7i2=7i;  tj2  = 


p^exp(ip-i) 

^Po 


l  +  r(z  -l)-k(z-l)/ z 
z  =  (l+fl  +  4ico  )/2;  i=V^;  pi=0.1pQ-,{#)  5=v-r-pk\ 

g,  =  {g2  +  [z^(l - £) - 1  +  ^(1  +  £)]}  +  v[C/, (2z2  - 2z - 1)  +  t/jz - Zi  - 1];  D  =  k^r-\-, 


gi 


_(^2 — - - — ')  +  (z2-Zj)(D - - — );  Zj  =(l+fl-4i(i)  yi',  Z2  -  (1+ Vl  +  )/2; 


26? 


_  .  v^here  E  is  activation  energy  of  gasification,  and  is  burning  surface  temperature  under 

/  E 


stable  combustion  at  pressure  po.  Functions  U\  and  Uj  can  be  found  by  eq.(7.4),  in  which  function  z  must 
be  substituted  by  Z\  and  Z2,  correspondingly. 

It  is  important  to  stress  the  new  aspect  of  the  problem:  the  complex  response  function  Vi  depends  on 
the  phase  shift  y/  between  the  first  and  the  second  harmonics.  Because  of  that  the  result  of  calculations 
was  received  for  four  values  oix]/  :  0, 7c/2,  jt,  and  37c/2. 


The  calculated  typical  nonlinear  amphtudes  Re{Vi}  of  burning  rate  response  functions  and  phase 
shift  of  the  response,  Im{Vi},  as  dependencies  on  a  are  shown  on  Figures  52-22  for  mixture  1.  Similar 

functions  were  obtained  for  other  mixtures.  Table  9  shows  values  of  Re{Vl}m  and  <»„  ( as - ^-^)  of 


the  functions.  Values  of  here  are  also  the  resonance  frequency,  corresponding  to  Re{Vl  }m. 

Relative  standard  deviations  of  values  Re{Vi}ni  obtaining  are  about  ±  (20-40%). 

It  can  be  seen,  that  nonlinear  correction  of  response  functions  strongly  depends  on  mixture  and  on 
phase  shift  yf.  At  ^=0:  values  of  Re{Vi}m  for  mixtures  1  and  4  are  close  to  linear  Re{U}n)  of  these 
response  functions  practically  at  all  studied  pressures.  The  values  of  mixture  3  a  bit  more  than  linear 
values,  but  mixture  2  has  a  very  large  increasing  of  Re{Vi}m  in  comparison  with  linear  Re{U}m  of  this 
mixture.  At  yf=  %ll\  values  of  Re{Vi}ra  for  all  mixtures  even  less,  than  values  of  linear  Re{U}m. 

At  y/=  7t:  values  of  Re{Vi}ni  always  more  than  linear  Re{U}ni  for  all  mixtures,  but  that  difference  for 
mixtures  1,  4  and  3  it  is  not  large.  However,  for  mixture  2  it  is  significant.  At  y/=  37t/2:  values  of 
Re{Vi}m  have  maximal  values  at  every  pressure.  However,  for  mixtures  1,  3  and  4  they  are  not  very 
large,  but  for  mixture  2  they  are  very  large:  for  example,  at  2  MPa  Re{Vi}m=3.38,  when  linear 
Re{U}n,=0.84. 
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Thus,  influence  of  y/  has  no  monotonic  character:  Re{Vi}ni  increases  with  phase  shift  at  %  and 
'inti,  but  at  7c/2  values  of  Re{Vi}m  are  small  for  all  mixtures,  and  even  less  than  values  of  linear 
Re{U}ni.  It  implies,  that  the  response  at  yr=  tUI  partly  compensates  responses  at  other  values  of  phase 
shift  (because  phase  shift  goes  through  interval  0-2%), 

Values  of  here  are  also  strongly  depending  on  y/.  As  a  rule,  values  of  On  at  y/=0  and  3%I2 

significant^  higher,  than  those  at  linear  approach,  but  at  y/=  %I2  and  %  the  values  about  two  times  less, 
than  at  linear  case. 

The  fiinal  conclusion  of  the  section  is  as  follows:  the  presented  data  for  values  of  Re{Vi}n,  show, 
that  the  obtained  for  linear  case  ‘‘sequence  of  preferences”  in  section  7. 1  is  valid  also  for  the  second  order 
interaction  of  the  first  and  second  modes. 

7.2.2.  The  second  order  self-interaction  of  the  first  harmonic. 

Obviously,  for  the  second  order  self-interaction  of  the  first  harmonic  we  have  the  same  equations  (7.5), 
(7.6).  After  presenting  of  the  expressions  for  dimensionless  pressure  and  burning  rate  in  complex  forms 
and  solving  of  the  no  stationary  heat-conduction  equation  in  the  condensed  phase,  the  following  solution 
can  be  obtained  for  linear  response  function  plus  nonlinear  correction  of  the  function  by  the  second  order 
self-interaction  of  the  first  mode: 


V2==U2+[/n  (7.8) 

where:  U2  -  a  linear  response  fiinction,  calculated  by  eq.(7.4),  in  which  z  was  substituted  by  Z2,  and 
U\\  1/1  ///2  is  the  member  of  nonlinear  correction,  i.e.  a  result  of  second  order  self-interaction  of  the  first 
harmonic.  Here: 


t/„  = 


Gj  +  G2  +  G3 


5=v.r-^.k, 


l  +  r(z-l)-A:(z-l)/z’  -  -  2p, 

z=(l+Vl  +  4/fi)  )/2;  i=V^  pi=0.1po,i#)  z^  =  (1+Vl  +  8/0 )/2;  D  =  k  +  r-l; 

+  C73  =v?/-[1  +  ^(z2-z)]; 


s  = 


G> 


Function  U can  be  calculated  by  eq.(7.4). 


The  response  function  (7.8)  also  depends  on  phase  shift  y/,  and  the  calculations  also  was  made  for 
four  values  of  yf  :  0,  jt/2,  ji,  and  3jt/2. 

The  calculated  typical  nonlinear  amplitudes  Re{V2}  of  burning  rate  response  functions  and  phase  of  the 
response,  Im{V2},  as  dependencies  on  0  are  shown  on  Figures  23-30  for  mixture  1.  Similar  functions 

were  obtained  for  other  mixtures.  Table  10  shows  values  of  Re{V2}m  and  (a„  (  as  )  of  the 

0) 

n 

functions.  Relative  standard  deviations  of  values  Re{V2}m  obtaining  are  about  ±  (20-40%). 

It  can  be  seen,  that  second  order  self-interaction  of  the  first  harmonic  gives  a  very  small  corrections  to 
linear  Re(U}ni  for  all  investigated  mixtures  and  for  all  phase  shifts  yi.  Only  eftect  of  decreasing  of  values 
0n  in  two  times  is  observed.  It  is  a  natural  consequence  of  self-interaction  of  the  first  mode. 

So,  the  analysis  of  this  section  also  show,  that  “sequence  of  preferences”  is  observed  also  for  the 
second  order  self-interaction  of  the  first  mode. 


7.2.3.  Increased  level  of  pulsating  pressures. 

The  calculations  of  previous  paragraphs  7.2.1  and  7.2.2  were  performed  at  small  amplitude  of  pressure 
pulsations;  p\  =  O.l  po.  It  is  important,  from  scientific  and  practical  points  of  view,  to  estimate  the 
influence  of  increased  level  of  pulsating  pressures  on  the  obtained  nonlinear  response  functions.  It  can  be 
made,  because  the  level  of  pressure  pulsations  is  included  in  the  calculations  in  evident  form;  -  see  above 
(#).  Of  course,  assumption  of  small  amplitudes  of//,  is  the  base  of  the  theory  of  nonlinear  response 
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functions,  and  because  of  that  the  estimations  at  increased  pulsating  pressures  may  give  only  a  tendency 
of  the  influence  on  the  response  functions.  In  this  section  we  are  estimating  this  tendency. 

Calculations  by  formulas  (7.7)  and  (7.8)  were  made  at  pulsating  pressure p\  =  0.3  po  and 
p\  =  0.6  po .  The  results  of  estimations  are  shown  on  Tables  9  and  10. 

It  can  be  seen,  that  influence  of  increased  p\  on  the  second  order  self-interaction  of  the  first  harmonic 
doesn’t  large.  Indeed,  the  increased  p\  gives  at  the  self-interactions,  as  a  rule,  about  10%  increase  of 
Re{V2}m  values,  and  doesn’t  change  the  character  of  the  response  fimctions.  Influence  of  pressure 
amplitude  p\  =  0.3  po  on  second  order  interaction  of  the  first  and  second  modes  is  more  significant. 
Indeed,  the  increased  p\  gives  here,  as  a  rule,  increased  values  of  Re{Vi}m  on  about  20-60%,  but  also 
without  principal  changing  of  character  of  functions  Re{Vi(t»)}  and  Im{  Vi(e))}. 

Thus,  we  can  conclude,  that  influence  of  elevated  oscillating  pressures  doesn’t  change  the  main  result  of 
section  7  about  “sequence  of  preferences”;  it  is  valid  also  for  elevated  pulsating  pressmes: 

[GAP/HNrW]<[(AMMO-BAMO)/HMX]<[GAP/HMX]<[(AMMO-BAMO)/HNIW]  -  NONLIN  (x) 

On  the  left  side  here,  as  above,  stands  the  minimal  response  and  on  the  right  side  -  the  maximal  one. 
Mark  ‘NONLIN’  show,  that  it  is  the  nonlinear  case  with  inclusion  of  the  case  of  elevated  oscillating 
pressures. 

So,  we  can  conclude  that  linear  approach  to  the  problem  of  response  function  of  the  investigated 
mixtures  is  a  good  approach,  and  nonlinear  approach,  including  elevated  pulsating  pressme  level,  gives 
only  not  large  numerical  corrections,  without  principal  changing  of  the  main  results. 


8.  CONCLUSIONS 

Main  Features  of  Combustion  Mechanisms  of  the  Studied  Mixtures 
At  Pressure  and  Sampie  Temperature  Variations 
and  Under  Conditions  of  Stabie  and  Osciiiating  Pressures. 

1)  A  significant  difference  in  burning  rates  of  mixtures  on  base  of  HNIW  and  on  base  of  HMX  takes 
place:  burning  rates  of  HNIW-based  mixtures  are  significantly  higher. 

2)  Influence  of  type  of  binder  on  burning  rate  of  HNIW-based  mixtures  is  significantly  higher  than  that  on 
burning  rate  of  HMX-based  mixtures.  (AMMO-BAMO)/HNIW  mixture  has  higher  burning  rate  than  that 
of  GAP/HNIW  mixture.  (AMMO-BAMO)/HNIW  mixture  has  also  especially  high  bimi-rate  temperature 
sensitivities  at  negative  sample  temperatures. 

3)  Temperatures  of  bimiing  surface  of  HNIW-based  mixtures  are  significantly  less  than  those  of  HMX- 
based  mixtures  at  all  investigated  pressures  and  sample  temperatures. 

4)  Both  type  of  mixtures  have  two-zone  structure  of  temperature  profiles  in  the  gas  phase  at  middle 
pressures. 

5)  Both  type  of  mixtures  have  wide  zone  reaction  in  the  gas  phase  and  distributed  (along  the  zone)  the 
functions  of  heat  release  rate  at  all  investigated  pressures  and  sample  temperatures. 

6)  Significant  heat  release  in  solid,  increasing  with  pressure,  and  relatively  small  heat  feedback  from  gas 
to  sohd  due  to  conduction  and  radiation,  decreasing  with  pressure,  are  observed  for  both  type  of  mktures. 

7)  Burning  rate  control  regions  in  the  combustion  waves  of  both  mixtures  are  placed  on  the  burning 
surface  and  in  the  gas  phase  close  to  surface,  at  low  temperatures;  high  temperature  gas  regions  cannot 
influence  the  burning  rate  at  all  investigated  pressures  and  sample  temperatures.  Both  mixtures  have 
stable  combustion  under  investigated  conditions.  Intensity  of  chemical  reactions  in  the  gas  phase  close  to 
the  burning  surfaee  for  HNIW-based  mixtures  is  significantly  higher  than  the  intensity  for  HMX-based 
mixtures. 

8)  The  investigated  mixtures  have  different  macrokinetic  laws  of  solid  gasifications,  depending  on 
oxidizer  and  practically  independent  on  binder.  Activation  energy  of  HNIW-based  mixture  gasification  is 
equal  to  32.8  Kcal/mole  and  that  of  HMX-based  mixtures  -  28  Kcal/mole  (at  Ts<500®C). 

9)  Pressure-driven  bum-rate  response  functions,  obtained  for  the  studied  mixtures  at  linear  and  nonlinear 
approximations  for  the  second  order  interactions,  have  complex  character  and  they  are  typical  for  many 
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condensed  systems.  The  obtained  results  allow  one  to  establish  the  best  and  the  worst  mixtures,  from 
point  of  view  of  the  reaction  of  burning  rate  on  pressure  pulsations;  the  best  mixture  has  the  smallest 
burning  rate  response  and  the  worst  mixture  -  the  largest  one.  The  reaction  was  established  by  values  of 
Re{V}„  -  of  the  real  part  of  the  response  functions  for  the  studied  mixtures.  On  the  base  of  the  values,  the 
following  «  sequence  of  preferences  »  was  suggested  for  the  studied  mixtures; 

[GAP/HNIW]  <  [(AMMO-BAMO)/HMX]  <  [GAP/HMX]  <  [(AMMO-BAMO)/HNIW] 

On  the  left  side  stands  the  minimal  response  and  on  the  right  side  -  the  maximal  one. 

The  «  sequence  of  preferences  » is  the  base  for  estimations  of  ‘combustion  quality’  of  new  mixture 
burning  in  field  of  pulsating  pressure. 

The  « sequence  of  preferences » obtained  by  the  linear  approach,  probable,  will  be  valid  also  for 
nonlinear  cases. 


Recommendations 

In  future  work  must  be  investigated  combustion  mechanisms  of  mixtures  on  base  of 
HNIW,  with  catalytic  additions,  correcting  burning  rates,  and  having  new  binders.  It  is  important 
to  perform  investigations  of  gas  velocity  driven  bum-rate  response  functions. 
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TABLES  OF  FINAL  REPORT  (item  OOOIAD) 


Table  1;  Burning  wave  parameters  of  mixture  1  (GAP/HNIW,  20:80, 1.77  g/cm^) 
at  pressure  p  and  sample  temperature  To  variations.  Presented  at  To,*C;  -100/20/100. 


N 

p,  MPa 

0.1 

0.5 

1.0 

2.0 

5.0 

8.0 

10.0 

1 

m,  g/cm^s 

-/0.28 

0.37/0.5 

0.58/0.8 

0.95/1.25 

1.8/2.05 

2.65/3.25 

3.3/4 

/0.35 

/65 

/0.95 

/1 .4 

/1.3 

13.65 

/~4.45 

Ts/C 

-/248 

256/270 

270/285 

290/303 

315/324 

330/344 

345/355 

■ 

/255 

/278 

1191 

/308 

/328 

/350 

1360 

B 

(p-10"'‘,  K/cm 

-/2/2 

3. 3/3/3 

4.5/4/3.5 

6.5/5.5/4.5 

10/9/7.5 

12.6/1 1/10 

14/13/12 

■1 

-/1 2 

15.6/10.5 

13.6/9 

12.3/8 

10.6/8 

9/6.6 

8.3/6.3 

4,  Gdi/g 

19.1 

/8.3 

/6.6 

,  /6 

fe.i 

/5.4 

/5.4 

5 

qr,  cal/g 

-ms.i 

5. 4/4/3 

2.5/2.572 

2.1/4/3.6 

7.8/6.8/6 

977.4/6.6 

7.3/877.2 

6 

Q,cal/g 

-/61/39 

104/73/51 

114/81/59 

122/87/63 

127/91/68 

133/99/76 

140/103/7 

7 

/,  pm 

80/70/55 

65/55/45 

48/45/40 

35/33/30 

30/27/20 

25/23/20 

8 

pm 

-/50/60 

40/40/40 

30/30/36 

20/25/32 

15/20/28 

12/17/20 

10/15/18 

9 

X■10^cmVs 

1. 7/2.0 

/2. 1/2.5 

2.6/3. 0 

3. 6/4.0 

4.2/5.0 

4.615.5 

/1 .4/1. 5 

/2.0 

11.5 

/3.2 

/4.0 

/4.5 

/5.0 

10 

Xs-10'', 

_ 

10/12 

13/16.7 

16/18.6 

23/24 

26/31 

28/34 

cal/cm-s-K 

78.7/9.3 

/12.5 

/1 5.5 

/20 

115 

/28 

731 

11 

Ti,“C 

-/800 

880/900 

980/1000 

1150/ 

-/71350 

-f-l- 

-/-/- 

/850 

/950 

/1 150 

/1250 

12 

H 

o 

O 

-/-/- 

-/2500 

2500/2600 

2630/2750 

2750/2850 

2850/2930 

2850/2930 

/2550 

11650 

/2800 

/2950 

/3000 

3000/ 

13 

Li,  mm 

-^>3/^5 

>5/2/1 .5 

2.5/1 .5/1. 3 

0.9/1. 0/1 .2 

-1-10.3 

-/-/- 

-/-/- 

14 

L,  mm 

././■ 

-14/2.5 

5/4/3 

3.5/3/4 

2.8/2.573.5 

1.8/2.0/3 

1.6/1 .8/2 

15 

^lm 

-/28 

22/18 

15/11 

10/7.8 

6.0/5. 2 

4.0/3.3 

3. 2/2.7 

/23 

/1 3 

/lO 

fl.O 

/5.0 

/3.0 

11.5 

16 

n 

-/107 

>220/111 

160/136 

90/128 

475/240 

450/303 

500/333 

>220 

/no 

7127 

/\61 

/43 

/500 

7750 

17 

Wo,  kcal/cm^s 

-/1. 9 

4.0/5.0 

8.7/1 1 

20.7/23 

60/63 

113/122 

156/177 

/2.3 

/6.6 

/1 1.2 

/21.4 

/58.7 

/125 

/1 83 

15 


Table  2;  Burning  wave  parameters  of  mixture  2  ({BAMO-AMMO}/HNIW,  20:80, 
1.76  g/cm^  at  pressure  p  and  sample  temperature  To  variations. 
Presented  at  T„,“C:  -100/-50  /20/100. 


N 

p,  MPa 

0.5 

1.0 

2.0 

5.0 

8.0 

10.0 

1 

m,  g/cm^s 

-10.1 

0.75/1.0 

1.15/1.55 

2.473.3 

3.6/46 

4.3/5.55 

/0.8/1.1 

71.271.55 

72.1/2.5 

/4.4/5.0 

76.076.9 

/6.9/~8 

2 

Ts,‘’C 

-/278 

280/290 

296/308 

326/340 

344/356 

3527365 

/285^98 

7302/312 

7325/332 

7360/365 

^75/382 

73827-390 

3 

K/cm 

-/4.4 

5.6/5 

7/6.5 

11/10.5 

i’'4.5/i3"7 . 

16.5/15.8 

/4/2.8 

/4.6/3.5 

76/4.6 

710/8 

713/12 

/1 5/1 4 

4 

q,cal/g 

-/ll 

13/9 

1177.8 

8.7/6 

7.875.8 

7. 7/5.7 

/9/4.6 

77/4.2 

75.4/3,5 

74.473.1 

74.3/3.5 

74.3/3.6 

5 

qr,  cai/g 

-/3 

2.7/2 

4/3.2 

5.8/4.2 

6.7/5.2 

6.5/5 

/2.5/1.8 

71.671.3 

72.4/2 

/3 .5/2.8 

7473.5 

/4.6/4 

6 

Q,  cal/g 

-/1 01 

117/108 

124/114 

135/126 

141/131 

144/135 

781/63 

790/69 

799776 

71 14/87 

7116/92 

7118/94 

7 

/,nm 

-750 

50/44 

40/36 

36/26 

22/21 

20/19 

745/48 

740/40 

734735 

725/24 

722/16 

720/17 

8 

-725 

15/22 

13/19 

10/14 

8/12 

12/10 

735/45 

730/33 

724/30 

722/20 

718/18 

/1 6/1 8 

9 

X-IO’,  cmVs 

-71.7 

2.172.5 

2.6/3 

4/4.8 

5/5.5 

5.6/6 

72/3 

72,7/3.5 

74/5 

75/5.5 

76/6 

/6.577 

m 

VI  o\ 

”7i'o . 

12/15 

16/19 

25/29 

28/33 

30/35 

cal/ctns-K 

712/16 

716/20 

722/28 

732/35 

735/40 

/40/43 

ii 

t,/c 

-7900 

850/950 

900/1000 

950/1050 

-/- 

-/- 

71000/1100 

71050/1150 

71100/1180 

/1 150/1200 

A/- 

/-/- 

12 

Tf,”C 

-72140 

2266/2300 

2500/2540 

2766/2770 

2800/2850 

2800/2850 

72200/2250 

72400/2500 

72600/2700 

72850/2950 

72930/3000 

/2930/3000 

i3 

Li,mm 

-A).8 

0.5/0.7 

6,5A).6 

0.2/0.3 

-A 

•:/- 

/1 .0/1 .5 

/0.9/1.2 

70.6/0.7 

/0.25/0.4 

AA 

/-/- 

14 

L,  mm 

15/3.5 

2/3.1 

ij8/2"8 . 

1. 5/2.1 

1.4/1 .7 

/5.0/5.0 

75.0/6.0 

/4. 5/5.0 

/4.0/4.5 

/3.0/3.5 

/2.0/2.5 

15 

-712 

11/9 

7.4/6 

. . 

3/2.4 

2.5/2 

712/9,4 

/8.4/6.9 

74.9/4.4 

llMl.l 

72.1/1.9 

/1 .8/1 .6 

16 

n 

-/67 

45/78 

68/100 

54/110 

566/875 

560/«50 

/83/160 

/1 10/1 70 

7122/160 

7100/180 

/1 430/790 

/1 100/940 

17 

-710 

14717 

27/34 

89/117 

177/214 

711/10 

719/18 

743/39 

7150/138 

/270/280 

360/380 

16 


Table  3;  Burning  wave  parameters  of  mixture  3  (GAP/HMX,  20:80, 1.74  g/cm^) 
at  pressure  p  and  sample  temperatnre  To  variations.  Presented  at  To:  -100"C/20"C/90"C 


p,  MPa 

0.5 

1.0 

2.0 

5.0 

10.0 

m,  g/cm^s 

-/O.IA 

0.155/0.21/0.26 

0.28/0.37/0.43 

0.64/0.77/0.88 

1.1/1.26/- 

T^“C 

-/330/- 

345/365/375 

380/400/410 

445/460/470 

485/492/- 

(p-lO-'.K/ctn 

-13.51- 

2.0/4.0/5.0 

3.5/5.0/7.8 

9.0/10/11 

11.5/12/- 

q,cal/g 

-mi- 

35/37/42 

29/31/43 

35/30/33 

27/25/- 

Or,  cal/g 

-10.11- 

0.8/0.7/0.6 

0.8/.06/0.6 

2.7/2.6/2.4 

3.1/3.71- 

Q,  cal/g 

-mi- 

152/120/90 

180/140/103 

215/164/140 

223/180/- 

/,  pm 

-lies/- 

150/108/106 

108/78/75 

50/58/54 

48/48/- 

/m,  pm 

-1451- 

30/40/55 

30/35/40 

40/30/38 

20/20/- 

X*10^cm^/s 

-n.oi- 

1.0/1 .3/1.6 

1.7/1 .7/1 .9 

2.3/2. 5/2.8 

3. 1/3.3/- 

cal/cm  s-K 

-16.11- 

6.7/9.0/1 1 

12/12/13 

16/17/19 

22/23/- 

Ti,  °C 

-/980/- 

-/1050/1000 

-f-inoo 

-/-/- 

-/-/- 

Tf/C 

-/1900/- 

1950/2050 

/2150 

2150/2200 

/2280 

2300/2420 

/2500 

2300/2420/- 

Li,  inm 

-16001- 

-/500/900 

-1-1350 

-f-f- 

-/-/- 

L,  mm 

-12.01- 

1.5/1 .7/2.3 

1.5/1 .6/2.3 

2.0/2.0/2.5 

1. 5/1.7/- 

-1431- 

43/24/19 

20/15/13 

21/8.0/7.0 

5. 6/4.9/- 

n 

-1151- 

35/96/120 

75/106/180 

95/250/360 

270/300/- 

Wo,kcal/cm^s 

-n.4/- 

1.3/3 .5/5.5 

4.1/7.8/14 

24/33/40 

54/64/- 

Table  4;  Burning  wave  parameters  of  mixture  4  ({BAMO-AMMO}/HMX,  20:80, 
1.72  g/cm*)  at  pressure  p  and  sample  temperature  T®  variations. 
Presented  at  T„:  -100”C/20"C/90'’C 


N 

p,  MPa 

0.5 

1.0 

2.0 

5.0 

10.0 

1 

m,  g/cm^s 

-/O.l/- 

0.16/0.220.28 

0.32/0.41/0.48 

0.72/0.85/0.98 

124/1.44/- 

2 

T,/C 

-/330/. 

345/365/380 

390/407/420 

460/470/480 

492/500/- 

3 

(plO-^,K/cm 

-/2.2/- 

3.4/3.7/5.0 

5.0/6.0/6.5 

8.0/9.0/11.0 

10.4/12.0/- 

4 

q,cal/g 

•/46/. 

44.6/37.0/41.0 

36.0/34.0/33.0 

27.0/26.0/29.0 

22.0/22.0/- 

5 

qr,cal/g 

-/0.6/- 

.06/0.6/0.4 

.05/0.4/0.4 

2.0/1 .8/1 .7 

3.4/4.0/- 

6 

Q,  cal/g 

-/96/- 

150/120/93 

180/140/120 

215/172/145 

230/185/- 

7 

/,nm 

-/195A 

142/110/100 

100/80/78 

45/58/55 

53/45/- 

8 

/m,nm 

-/35/- 

35/34/40 

30/35/40 

25/30/45 

20/30/- 

9 

X*10^  cm  Vs 

-/I.IA 

1.2/1 .4/1 .6 

1.7/1. 9/2.1 

2.7/2.8/3.1 

3.6/3.8/- 

10 

Xs*10\  cal/cmsK 

-/7.6/- 

8.3/9.6/11.0 

11.7/13/14 

18.5/19.3/21 

24.8/26/- 

11 

Ti/C 

-/900A 

-/1020/1050 

-/-/- 

-/-/- 

-/-/- 

12 

Tf/C 

-/1 850 
/- 

2000/2100 

2200 

2200/2300 

/2350 

2280/2380 

/2450 

2300/2400 

/- 

13 

Li,  mm 

-/500/- 

V400A700 

-/-/- 

-/./- 

V-/- 

14 

L,  mm 

./2.2/- 

1.9/2.0/2.3 

2.0/2.2/2.3 

1.7/2.0/2.2 

l.l/l.O/- 

15 

d,  ^im 

-/51A 

32/24/20 

17/13/12 

8.3/7.0/6.2 

4.9/4.0/- 

16 

Q 

-/9.8/. 

53/17/35 

100/154/167/ 

145/258/323/ 

204/225/- 

17 

Wo,  kcal/cm^s 

-/0.9/- 

2.2/3 .4/6.0 

6.7/10/13 

24/33/46 

55/73/- 

Table  5i  Coefficients  of  heat  conduction  Xg  and  specific  heat  Cp  of  the  gas  phase 


mixtures 

T/C 

300 

500 

800 

1000 

1500 

GAP/HNIW 

X,g*10'’,  cal/cm  c-K 

1.6 

2.1 

2.8 

3.3 

4.4 

Cp,  cal/g-K 

0.33 

0.34 

0.36 

0.37 

0.40 

(BAMO- 

AMMO)/HNIW 

Xg’lO^,  cal/cm  C'K 

1.6 

2.1 

2.8 

3.4 

4.5 

Cp,  cal/g  K 

0.33 

0.34 

0.36 

0.37 

0.40 

GAP/HMX 

Xg‘10^,  cal/cm  c  K 

2.0 

2.6 

3.4 

4.0 

5.6 

Cp,  cal/g  K 

0.41 

0.43 

0.45 

0.46 

0.50 

(BAMO- 

AMMO)/HMX 

Xg-10^,  cal/cm-c-K 

1.9 

2.5 

3.3 

3.9 

5.5 

Cp,  cal/g  K 

0.41 

0.43 

0.45 

0.46 

0.50 

Table  6:  Sensitivities  of  burning  rate  and  surface  temperature 
to  pressure  and  sample  temperature  of  the  studied  mixtures. 
Sample  temperature  To=20"C 


Systems 

P,  IVIPa 

0.1 

0.5 

1 

2 

5 

8 

10 

Activation 
energy  of 
gasification 

Mixture  1 

GAP/eNIW, 

property 

0.28 

0.63 

0.37 

0.088 

0.06 

521 

0.31 

0.78 

0.50 

0.106 

0.076 

543 

0.23 

0.60 

0.60 

0.102 

0.09 

558 

0.19 

0.54 

0.60 

0.08 

0.084 

576 

0.13 

0,40 

0,80 

0,06 

0,08 

597 

0.17 

0,55 

0,93 

0,09 

0,14 

617 

0.15 

0,50 

0,93 

0,075 

0,136 

628 

p,%/K 

k 

V 

r 

Ts,K 

£,  Kcal/mol 

32.8 

P,%/K 

0.28 

0.30 

0.36 

0.30 

0.30 

0.28 

Mixture  2 

k 

- 

0.74 

0.85 

1.10 

1.02 

1.12 

1.01 

V 

- 

0.55 

0.70 

0.82 

0.73 

0.68 

0.68 

(AMMO- 

r 

- 

0.13 

0.15 

0.15 

0.18 

0.18 

0.20 

32.8 

BAMO)/BDVIW, 

i» 

0.10 

0.12 

0.115 

0.11 

0.10 

0.10 

Ts,K 

- 

558 

575 

598 

633 

648 

655 

Mixture  3 

p,%/K 

0.28 

0.23 

0.18 

0.17 

0.17 

k 

- 

- 

0,97 

0,87 

0,79 

0,787 

0,80 

GAP/BMX, 

v 

- 

- 

0,83 

0,80 

0,76 

0,76 

0,76 

r 

- 

- 

0,17 

0,155 

0,13 

0,10 

0,09 

28.0 

_ 

- 

0,16 

0,165 

0,115 

0,104 

0,092 

Ts,K 

- 

- 

638 

673 

733 

756 

765 

Mixture  4 

p,%/K 

“ 

0.28 

0.20 

0.16 

0.13 

0.12 

k 

- 

- 

0,95 

0,759 

0,72 

0,61 

0,576 

(AMMO- 

v 

- 

0,90 

0,82 

0,80 

0,79 

0,79 

BAMO)/HMX, 

r 

- 

- 

0,19 

0,16 

0,105 

0,066 

0,067 

28.0 

20:80 

P 

- 

- 

0,174 

0,168 

0,11 

0,094 

0,092 

Ts,K 

.  - 

- 

638 

680 

743 

763 

773 

Table  7:  Dependencies  of  sensitivities  for  burning  rate  and  surface  temperature 
of  system  (AMMO-BAMO)/HNIW,  20:80,  at  negative  sample  temperatures. 

Criteria  of  stable  combustion. 


0.5 

1 

2 

5 

8 

10 

-100 

-50 

-100 

-50 

-100  1 

-50 

-100 

-50 

0.20 

0.56 

0.44 

0.47 

0.47 

0.42 

0.42 

0.42 

0.42 

0.40 

0.40 

0.66 

2.13 

1.5 

1.87 

1.68 

1.80 

1.64 

1.80 

1.70 

1.37 

1.25 

0.10 

0.20 

0.19 

0.20 

0.23 

0.22 

0.25 

0.24 

0.25 

0.27 

0.25 

2.0 

1.32 

0.7 

1.03 

0.55 

0.93 

0.5 

0.2 

0.1 

19 


Table  8:  Response  function  characteristics  of  tfie  studied  mixtures. 
Linear  approximation,  T«=20®C. 


Properties 


Systems 

P> 

MPa 

Re{U}„ 

e>n 

KHz 

Re{U}20 

Im{U} 

niax+ 

e>im+ 

®^in+5 

KHz 

lm{U} 

min 

^im 

KHz 

0.1 

0.55 

6.75 

0.12 

0.46 

0.07 

1.0 

0.018 

-0.12 

15 

0.28 

C'  Al>/ 

mn 

0.86 

5.75 

0.23 

0.70 

0.14 

1.25 

0.05 

-0.24 

16 

0.63 

XrAr/ 

HNIW 

1 

5.5 

0.45 

0.73 

0.11 

1.0 

0.08 

-0.17 

17 

1.38 

2 

0.84 

6.0 

0.98 

0.76 

0.10 

1.0 

0.16 

-0.14 

18 

2.93 

5 

1.04 

7.75 

2.56 

1.00 

0.10 

1.0 

0.33 

-0.08 

18 

5.94 

8 

1.30 

6.0 

4.02 

1.18 

0.15 

1.0 

0.67 

-0.24 

18 

12.1 

10 

1.27 

6.0 

5.52 

1.20 

0.14 

1.0 

0.92 

-0.20 

18 

16.6 

(AMMO- 

Ti  k’XMCWI 

4.5 

0.45 

0.13 

1.0 

0.10 

-0.26 

17 

1.72 

1 

1.21 

4.75 

0.81 

0.80 

0.21 

1.25 

0.21 

-0.39 

18 

3.1 

BAMUj/ 

HNIW 

2 

2.03 

6.0 

2.12 

1.20 

0.49 

1.75 

0.62 

-0.80 

17 

6.0 

5 

1.51 

4.75 

5.9 

0.80 

0.32 

1.5 

1.88 

-0.52 

15 

18.8 

8 

1.62 

5.0 

9.7 

0.80 

0.39 

1.75 

3.4 

-0.60 

14 

27 

10 

■B 

10.6 

0.80 

0.29 

1.25 

2.95 

-0.40 

13 

30 

1 

H 

0.05 

0.90 

0.29 

1.25 

0.014 

-0.60 

15 

0.17 

GAP/ 

HMY 

2 

1.36 

4.5 

0.12 

0.85 

0.22 

1.0 

0.026 

-0.50 

15 

0.39 

5 

1.28 

5.25 

0.40 

0.90 

0.22 

1.25 

0.096 

-0.36 

16 

1.23 

JXiVLA 

8 

1.34 

6.5 

0.83 

0.93 

0.23 

1.5 

0.19 

-0.34 

16 

2.1 

10 

1.40 

7.5 

1.18 

1.20 

0.26 

1.5 

0.24 

-0.32 

18 

2.8 

(AMMO- 

1 

1.62 

4.0 

0.05 

0.95 

0.29 

1.0 

0.014 

-0.59 

15 

0.20 

2 

1.26 

4.0 

0.12 

0.95 

0.18 

1.0 

0.03 

-0.40 

15 

0.45 

HMX 

5 

1.31 

6.0 

0.52 

0.80 

0.21 

1.25 

0.11 

-0.30 

18 

1.55 

8 

1.22 

8.0 

1.18 

0.79 

0.17 

1.25 

0.19 

-0.18 

18 

2.66 

10 

1.18 

7.75 

1.44 

0.79 

0.16 

1.0 

0.19 

-0.16 

20 

3.72 

20 


Table  9.  Response  function  characteristics  at  nonlinear  approximations: 
Second  order  interaction  of  the  first  and  second  harmonics. 
Values  of  Re{Vi}m/  are  presented.  To=20®C. 

Parameter  n  is  a  coefficient  in  expression:  pi=n^p0  ; 


p> 

MPa 

0.1 

0.5 

1.0 

2.0 

5.0 

10  1 

n 

0.1 

0.3 

0.1 

0.3 

0.1 

0.3 

0.1 

0.3 

0.1 

0.3 

0.1 

0.3 

Systems 

¥ 

0 

0.53 

0.70 

0.87 

1.1/ 

0.85 

1.0 

0.81 

0.90 

1.0 

0.96 

1.23 

1.40 

/8.0 

/20 

711 

20 

/II 

/20 

/1 4 

/20 

/20 

/20 

120 

/20 

n/2 

0.51 

0.46 

0.73 

0.66 

0.78 

0.73 

0.77 

0.72 

0.96 

0.93 

1.1 

1.1 

GAP/ 

/4.0 

/1.0 

/3.0 

/2.0 

/3.0 

/1.0 

/3.0 

/1.0 

/3.0 

/1.0 

/2.0 

/1.0 

HNIW 

7t 

. 

0.90 

... 

0.94 

1.1 

0.90 

_ 

1.15 

1.38 

1.4 

1.8 

/4.0 

/4.0 

/2.0 

/5 

/5.0 

/4.0 

/5.0 

/3.0 

3n/2 

. 

1.0 

_ 

1.0 

1.3 

0.97 

_ 

1.22 

1.64 

1.55 

2.2 

/8.0 

/8.0 

/10 

/10 

/20 

/20 

/1 2 

/1 8 

0.86 

1.1 

1.3 

1.76 

2.9 

5.1 

1.8 

2.7 

1.55 

2.2 

u 

/8.0 

/1 5 

/9.0 

/1 6 

/10 

/II 

/8.0 

/10 

18.0 

/10 

tz/2 

0.70 

0.70 

0.98 

0.92 

1.33 

1.34 

1.1 

1.0 

1.0 

0.94 

(AMMO- 

/3.0 

/1.0 

/2.0 

/1.0 

/2.0 

/1.0 

/2.0 

/1.0 

12.0 

/1.0 

BAMO)/ 

n 

. 

_ 

0.92 

_ 

1.32 

1.63 

2.4 

_ 

1.7 

2.1 

1.46 

1.83 

HNIW 

/3.0 

/3.0 

/3.0 

/4.0 

/3.0 

/3.0 

/3.0 

12.0 

3n/2 

. 

1.0 

_ 

1.54 

2.23 

3.4 

_ 

2.1 

3.3 

1.82 

2.75 

/6.0 

/6.0 

/6.0 

/6.0 

/5.0 

/5.0 

/5.0 

/5.0 

n 

1.8 

2.8 

1.5 

2.14 

1.36 

1.9 

1.55 

2.3 

u 

/8.0 

/10 

/8.0 

/1 2 

/1 2 

/1 9 

/1 7 

/20 

Tt/2 

. 

_ 

1.16 

1.13 

1.1 

1.0 

1.05 

1.0 

1.1 

1.1 

GAP/ 

n.o 

/1.0 

/2.0 

/1.0 

/2.0 

/1.0 

/3.0 

/1.0 

HMX 

K 

. 

_ 

_ 

1.77 

_ 

1.5 

1.9 

1.42 

- 

1.5 

2.1 

/3.0 

/3.0 

/2.0 

/4.0 

/5.0 

/4.0 

3k/2 

- 

_ 

_ 

2.23 

' 

1.8 

2.7 

1.65 

- 

1.9 

3.0 

/5.0 

/5.0 

/6.0 

/7.0 

/10 

/II 

n 

1.87 

2.87 

1.3 

1.74 

1.36 

1.9 

1.16 

1.3 

u 

/8.0 

/10 

/9.0 

/1 3 

/1 6 

/20 

120 

120 

(AMMO- 

n/2 

- 

_ 

_ 

- 

1.2 

1.17 

1.0 

0.95 

1.1 

1.05 

1.0 

1.0 

BAMO)/ 

n.Q 

/1.0 

/2.0 

/1.0 

12.0 

/1.0 

/3.0 

/1.0 

HMX 

n 

. 

_ 

1.8 

_ 

1.4 

1.7 

1.46 

_ 

1.32 

1.64 

/3.0 

/3.0 

I2.Q 

/4.0 

/5.0 

/5.0 

3n/2 

. 

_ 

_ 

2.23 

_ 

1.6 

2.22 

1.67 

- 

1.46 

2.14 

/5.0 

/5.0 

/6.0 

/8.0 

/1 4 

120 
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Table  10.  Response  function  characteristics  at  nonlinear  approximations: 
Second  order  self-interaction  of  the  first  harmonic. 

Values  of  Re{Va}m/  co„are  presented.  To=20®C. 

Parameter  n  is  a  coefficient  in  expression:  pi=n'Po ; 


11 

0.5 

1.0 

2.0 

5.0 

10 

n 

0.1 

0.6 

0.1 

0.6 

0.1 

0.6 

0.1 

0.6 

0.1 

0.6 

0.1 

0.6 

Systems 

¥ 

0 

0.51 

0.53 

0.86 

0.90/ 

0.88 

0.95 

0.86 

0.92 

1.1 

1.17 

1.9 

1.44 

12.0 

/2.0 

/3.0 

3.5 

/3 

/3.25 

/3.5 

/3.5 

/4.0 

/5.0 

72.75 

/4.0 

11/2 

0.51 

0.52 

0.86 

0.89 

0.88 

0.90 

0.85 

0.87 

1.1 

1.1 

1.27 

1.32 

GAP/ 

12.0 

12.0 

/3.0 

/3.25 

/2.75 

/3.0 

/3.0 

/3.25 

/4.0 

/3.75 

/3.0 

7Z.0 

HNIW 

n 

- 

- 

0.85 

- 

0.86 

- 

0.83 

1.02 

. 

1.24 

/2.75 

/2.5 

/3.0 

/4.0 

12.75 

371/2 

- 

- 

0.85 

- 

0.86 

0.84 

1.04 

1.26 

/3.0 

/2.5 

/3.0 

/4.0 

/3.0 

0 

_ 

0.87 

0.92 

1.21 

1.3 

2.0 

2.0 

1.52 

1.57 

1.38 

1.43 

/2.0 

/3.0 

12.0 

12.5 

/3.0 

/3;9 

12.5 

72.7A 

7225 

7215 

n/2 

- 

- 

0.86 

0.90 

1.22 

1.3 

2.05 

2.3 

1.53 

1.67 

1.38 

1.5 

(AMMO- 

12.0 

/3.0 

/2.5 

12.5 

/3.0 

/5.0 

12.5 

/3.0 

72.25 

72.75 

BAMO)/ 

n 

- 

- 

0.85 

. 

1.2 

2.0 

1.5 

1.35 

HNIW 

/2.5 

/2.25 

/3.0 

72.25 

72.25 

3n/2 

- 

- 

0.85 

- 

1.2 

_ 

2.0 

1.48 

1.34 

/2.25 

/2.25 

/3.0 

12.25 

72.25 

0 

- 

- 

- 

_ 

1.57 

1.66 

1.37 

1.47 

1.30 

1.40 

1.42 

1.49 

/2.25 

/2.75 

12.25 

/2.75 

72.75 

/3.25 

/3.75 

/4.5 

n/2 

- 

- 

- 

- 

1.58 

1.73 

1.37 

1.48 

1.30 

1.37 

1.42 

1.50 

GAP/ 

/2.25 

12.75 

/2.25 

12.5 

72.75 

/3.0 

/4.0 

74.0 

HMX 

n 

- 

- 

- 

- 

1.54 

- 

1.34 

1.26 

1.40 

12.25 

/2.0 

/2.5 

7Z.5 

3n/2 

- 

- 

- 

- 

1.53 

- 

1.34 

1.27 

1.40 

/2.25 

12.0 

12.5 

/3.75 

0 

- 

_ 

1.64 

1.77 

1.28 

1.40 

1.33 

1.44 

1.20 

1.31 

12.25 

/3.0 

/2.25 

12.75 

/3.0 

/3.75 

/4.0 

74.75 

(AMMO- 

n/2 

- 

- 

- 

- 

1.65 

1.82 

1.0 

1.35 

1.32 

1.39 

1.19 

1.23 

BAMO)/ 

12.0 

12.5 

12 

12.25 

/3.0 

/3.5 

/4.0 

74.0 

HMX 

n 

- 

- 

- 

- 

1.6 

. 

1.4 

- 

1.28 

1.16 

12.25 

/3 

/3.0 

/3.5 

3n/2 

- 

- 

- 

- 

1.6 

1.6 

_ 

1.29 

1.176 

12.0 

15 

/3.0 

7Z.5 
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FIGURES  OF  FINAL  REPORT  (item  OOOIAD) 


distance  x,  mm 


Figure  1:  Averaged  temperature  profiles  T(x)  of  mixture 
GAP/HNIW,  20/80,  at  0.1  MPa  and  different  To. 


23 


temperature  TC  temperature  T'C 


temperature  T'C 


25 
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Figure  7:  Averaged  temperature  profiles  T(x)  of  mixture 
GAP/HNIW,  20/80,  at  10  MPa  and  different  To. 


dimensionless  frequency 

Figure  8:  Linear  response  function  of  GAP/HNIW 
at  0.1  MPa.  To=20'C. 
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Figure  13:  Linear  response  function  of  GAP/HNIW 
at8MPa.To=20'C. 
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dimensionless  frequency 


Figure  14:  Linear  response  function  of  GAP/HNIW 
at  10  MPa.  To=20'C. 
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nonlinear  response  functions  Vi  nonlinear  response  functions  Vi 


Figure  20:  Nonlinear  response  functions  Vi  of  GAP/HNIW,  at  phase 
shift  i|i  =Tr  and  ip  =31x72;  p=5  MPa,To=20’C 
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Figure  22:  Nonlinear  response  functions  Vi  of  GAP/HNIW,  at  phase 
shift  (u  =TT  and  iii  =3Tr/2;  p=10  MPa,To=20'C 
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nonlinear  response  functions  V2 
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Figure  25:  Nonlinear  response  functions  Vz  (linear  functions  plus 
corrections  by  second  order  self-interaction  of  the  first  mode)  of 
GAP/HNIW,  at  phase  shift  ip  =0  and  qj  =tt/2;  p=1  MPa,To=20'C. 
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Figure  26:  Nonlinear  response  functions  Vz  of  GAP/HNIW,  at  phase 
shift  i|j  =TT  and  qi  =3tt/2;  p=1  MPa,To=20'C. 
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Figure  27:  Nonlinear  response  functions  V2  (linear  functions  plus 
corrections  by  second  order  self  interaction  of  the  first  mode)  of 
GAP/HNIW,  at  phase  shift  ip  =0  and  ip  =n/2:  p=5  MPa,To=20'C. 
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Figure  28:  Nonlinear  response  functions  of  GAP/HNIW 
ip  =TT  and  ip  =3Tr/2:  p=5  IVPa,To=20’C. 
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nonlinear  response  functions  V2 
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Figure  29:  Nonlinear  response  functions  Mi  (linear  Hinctions  plus 
corrections  by  second  order  self-interaction  of  the  first  mode)  of  GAP/HNIW, 
at  phase  shift  ip  =0  and  ip  =tt/2;  p=10  MPa,To=20'C. 
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Figure  30:  Nonlinear  response  functions  V2  of  GAP/HNIW,  at  phase  shift 
i|i  =TT  and  i|i  =3Tr/2;  p=10  MPa,To=20'C. 
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